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My aim will be threefold: 1) for each stage of stellar 

evolution, I will display the formula for the energy generation 
rate*; 2) I will discuss the uncertainties attached to each 
formula — uncertainties issuing both from the nuclear physics 
aspects of the problem and from possible variations in the 
chemical composition of the reacting material; 3) I will try to 
evaluate the size of the errors introduced in the analysis of 
stellar evolution when some minor nuclear reactions are 
neglected. (In other words, I shall try to find when these 
reactions can and cannot be neglected.) Because nuclear 
reactions produce new elements which may themselves become fuel 
at higher temperature, I shall consider some aspects of nucleo- 
synthesis. 

The basic equation describing the processes of energy 
generation and dissipation is 
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♦Much of the material presented here is discussed in Reeves (1964) 
where references are given. Only new material is given in any 
detail. 



This equation shows that the photon luminosity (per unit gram) , 

dL /dM , can be created by nuclear energy processes e M , by 
r r LN 

doing work on the gas or by decreasing its internal energy u. 

The term e v expresses the rate of energy dissipation by neutrinos. 
It may be worth mentioning that this formula implies the four 
basic interactions of modern physics: electromagnetic, nuclear, 

weak (neutrino emitting), and gravitational. In this talk I 
will consider only the nuclear and weak-interaction aspect of 
stellar energy generation. 

In equation (1) I have grouped together the terms e N and 
e . This way (e N - e v ) becomes a kind of "effective” nuclear 
energy rate and can be treated as such along the lines of 
previous analysis of stellar interiors (e.g., Schwa rzschi Id, 1957) . 

This term may become locally negative in a given star. 

During a nuclear burning stage this term must be mostly positive, 
or, more exactly 

J (e - e v ) dM > 0 , (2) 

M 

as we cannot otherwise expect the gravitational energy term to 
vanish, hence the contraction to stop. 

In cases where (e N - e v ) actually does become negative the 
luminosity will assume the following shape 
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In the zone A -* B photons are effectively being converted into 
neutrinos . 

To determine the presence of such zones it is most useful 
to work with two parameters 

n = d-tn e / d-tn T, (3) 

m = d-tn e / d-tn p . (4) 

n is called the temperature exponent and m the density exponent . 
Indeed we can write 

e (T, p ) = e 0 (T 0 ,p 0 ) (T/T 0 ) n (p/p 0 ) m (5) 

where T 0 and p Q are chosen in the range of interest for a given 
model. In this work I will quote the values of n and m for 
the various processes involved. 

The behaviour of the late stages of stellar evolution is 
ultimately related with the possibility (or impossibility) of 
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creating neutrinos in process involving a direct coupling 
between electrons and neutrinos, e.g., in processes like 
e” + y -* e” + "v or e + + e -* v + V . The physicists at 
CERN are now (1963) performing an experiment which may show 
the existence of the VT~ particle: a charged vector boson 

expected to mediate the weak interactions. (I use the terms 
"may show" as some experts in the field have emphasized the 
difficulty of drawing any definite conclusions from such in- 
volved and complicated experiments.) The existence of the W~ 
would most likely imply the reality of these neutrino processes. 
As of today however the subject is still unsettled. 

The presence (or the absence) of this extra mode of energy 
dissipation would alter considerably the time scales associated 
with the late stages of stellar evolution. Time scales may be 
evaluated (e.g., from star counts in various regions of the 
H.R. diagram) . Hence observational results may eventually be 
used as critical tests for (or against) the reality of the 
neutrino processes under consideration. 

In the present text, following the more or less general 
consensus of opinion, we shall assume the existence of these 
processes. Occasionally however we shall quote results which 
were computed without them, and discuss the difference. 



Three types of neutrino producing mechanisms are expected 
to be important during the course of stellar evolution: plasma 
neutrinos, photoneutrinos and pair-annihilation neutrinos*. 
Their respective domains are mapped in Figure (1) together 
with the profile of iso- intensity. We shall discuss each of 
them in turn later. 

The effect of these processes, together with the re- 
quirement stated in eqn. (2) seems to limit to four the 
maximum number of stellar nuclear-burning stages: hydrogen, 

helium, carbon, and oxygen burning. The time scale for neon 
photodisintegration is so short that the energy generated will 
at best slow down the contraction. (Nuclear reactions 
posterior to the fusion of oxygen will never be large enough 
to balance the tremendous output of neutrinos. 
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*Sakashita and Nishida (to be published) have called attention 
to another process: neutrino pair emission from excited nuclei. 

However its emission rate does in no condition appear to become 
as high as the sum of the rates of the three processes men- 
tioned here. 


* 


* A 

* V 


Through the work of Hayashi (1962) and Salpeter (to be 
published) we have learned that all stars will not pass through 
the four stages. Stars with masses smaller than about the 
Chandrasekhar limit (when due corrections are applied) will 
give up somewhere along the path to become dwarfs. The 
limiting lower masses are about 0.1, 0.5 and 0.7 to qualify , 

for Hydrogen, Helium and Carbon burning stages respectively. 

L 


HYDROGEN- BURNING- REACT IONS 


Here we consider a number of reactions competent in 
burning four protons into one Helium^ nucleus; yielding 6.68 
MeV per nucleon (1 MeV/nucleon ~ 10^® erg/gm) . The neutrino 
term of eqn. (1) is, at best a few percent of e^, so that 
(e^ - e^) never becomes negative. f 
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Proton- Proton Cycle 


The energy generation rate is given by 

e PP - e < x 4 = °> *,*<“."> • ( 6 ) 

Here e (X 4 = 0) is the rate of energy production if the core is 

devoid of Helium 4 , (X 4 = 0). (X ± will always be the fractional 

mass abundance of an isotope of mass number i) . ¥ then 

PP 

represent the effect of reactions induced by the He 4 nuclei 
themselves . 

Numerically we have 

e ( X 4 = 0) = (2.06 1 0.2) x lO 1 ^ f g pX 2 

1 i 1 1 , 1 1 

exp (- 33. 810 T,“ ^ 3 ) T~, 2 ^ 3 ( 7 ) 

b 6 

= T in millions degrees °K. 

The factor f^ ^ is the electron screening factor. For stars 

with M > 0.6 M Q a very good approximation is f.^ ^ = [l + 0.25 
1/2 3/2 

P /Tg ] . For smaller stars more involved formalism must 
be used. 

The factor g^ ^ is another correction term whose value 
is, close to unity. It grows from the value 1.02 at T c = 1, to 
1.09 at Tg = 15, to 1.19 at T^ = 50. The factor Ypp(a,w) can 
be written down as 

V p|p (a ' W) = { 1+Y [0.96 - 0.49 (W/l + W) ] } (8) 
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_ _ < 2; for T, > 30, 1.46, 

p» p b p, p 


We have 1 < Y 

The expression for y is 


j i 
< 


y = [ (1 + 2/a) ' - 1] a 


1/2 


(9) 


The value of y grows from zero to one as the temperature in- 
creases. Then a, the term directly dependent on He 4 abundance, 
has the form 


a = 5.48 x 10 


17 


(X /4X ) 
4 1 


exp (- 


100 T 1/3 ) 
6 


and finally, 


( 10 ) 


W = 1.22 x 10 16 (1 + p 1/2 /’I c 3/ ' 2 ) T~ 1/6 (1 + jLr 1 

6 o X x 

exp(- 102.6 T -1 ^) 

6 


( 11 ) 


The density exponent m = 1, and the temperature exponent are 

1/3 

given with quite good accuracy hy n = (11.3/Tg - 2/3). At 

T^ = 15, n = 4 while at T^ = 30, n ~ 3. In cases where the 
electron screening becomes important (small masses) the expo- 
nents can be altered to an appreciable extent and the above 
formula is no longer valid. 

The governing reaction in the proton-proton cycle is the 
reaction H' + H' -*D + e + + v. This reaction can be thought 
of as going into two steps: a) formation of an excited 
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L 
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diproton; b) beta decay of the diproton into a deuteron. The 
main uncertainties come from the second part, as both the beta 
decay coupling constant and the matrix element for the transi- 
tion are still slightly in doubt. The uncertainty or the rate 
is about 10%. 

We must now consider a complex of minor reactions which 

are involved in the proton-proton cycle. First there is the 
3 3 

He + He . Our estimate is based from an experiment done at 
Oak Ridge in 1954. The people involved have assigned a 50% 
accuracy to their results. However the people at the Calif- 
ornia Institute of Technology have their doubts about that. 

They plan to redo the experiments. Fortunately because of the 
form in which this cross section comes in the formalism of the 
energy generation, even a large change would not alter the 
total rate to an appreciable extent. In the absence of a 

better estimate I will use the Oak Ridge value and uncertainty. 

3 4 

The rate of He + He , through a recent experiment of 

Parker and Kavanagh (1963) is known to 15%. And the rate of 

the Be + p experiment is known to about 50%. These are the 

main uncertainties to worry about in the proton-proton cycle. 

These effects on the total rate are rougly as follows; 

For 1 <■ T < 10 only the p + p reaction matters (~ 10%) . In 
6 

the region from 10 < T < 20 and most sensitively around 

6 



3 o 3 j 

T ~ 14 the reactions He + He J and He + He* bring in an extra 

6 — 

source of uncertainty. The term a is known to about 60%. The 
term y is relatively insensitive to as in the worst case 
Ay/y ~ .20 (at a ~ 3) . In this region W « 1 so that 
y ~ 1 + Y is known to better than 10%. Hence in the range 

r / r 

10 < Tg < 20 the rate is known to better than 15%. 

In the next region 20 < Tg < 30, we have to consider the 

effect of the Be^ + p reaction. The expression W in eqn. (11) 

is known to about 50%. Its influence is most sensitively felt 

around w ~ 1 (T,. ~ 24) . It brings on Y and uncertainty of 

— b — pp 

about 4% (the effect of a is nil in that region) . Here again 

the rate is known to better than 15%. 

As you move toward higher temperatures Y becomes 

PP 

constant, and only the p + p uncertainty matters. 

Errors in the evaluation of the chemical composition also 

influence our knowledge of the rates. The proton-proton 

2 

reaction is proportional to (the hydrogen composition) but is 

also indirectly related to (the helium composition) through 
the term a in eqn. (9). If you change X^ by a factor of two 
you never change Ypp by more than 10 %. 

Carbon Nitrogen-Oxygen Cycle 


Here the energy generation formula is 


10 


( 12 ) 


e CNO " < 7 ' 9 - °- 8 > * 1027 f 14,l *14, 1 PX 14 X 1 V 2/3 


- 1/3 

exp(- 152.31 Tg ) erg/gra/sec. 

with f - [1 + 1.75 p 1/2 /T c 3/2 ], g = Cl + 10" 3 (3T 1/3 
14 , 1 6 14 , 1 o 

-4T 2 ' 3 )] . X, . is the fractional mass abundance of N, .. 

6 14 14 

The density exponent remains small (m = 1) while the 

1/3 

temperature exponent (n ~ 50.8/T - 2/3) becomes larger than 

O 

ten in all cases of interest (for the sun n ~ 19.5). For stars 
which draw most of their energy from this cycle such a high 
exponent drives the core into a convective state of energy 
transport. Such stars have central temperatures > 18. 

(For the sun Tg = 16 and the energy contribution of this cycle 
is only 3%. ) 

At temperatures T > 16 the cycle essentially transforms 

D 


all the isotopes of carbon and nitrogen into N 


14 


Hence the 


i 


term in eqn. (12) is very closely equal to the initial 

abundance of these isotopes (Caughlan and Fowler, 1963) . The 
same authors also show that at higher temperature the minor 
reactions 0^(p,y) F 17 (e + v) 0 17 , 0 17 (p,a) N 14 will effectively 

*1 fl 

throw in the game the 0 XO isotopes, thereby increasing corres- 
pondingly the value of X . Large uncertainties in the 
C> (p,a) N 1 rate make it difficult to evaluate with precision 
the onset of the transition. 
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As one approaches the very end of hydrogen-burning, a 
large number of minor reactions have to be taken into account. 
They have been considered by Parker et al (1964) . Although 
their energetic contribution is usually very small, it may 
become crucial in short periods. Great care must be taken in 
handling these ends of stages. 

After the exhaustion of hydrogen in the core, the nuclear 
burning of hydrogen takes place in a shell surrounding the 
core. The densities and temperature relevant to hydrogen 
burning both in the core and in the shell are pictured in 
Fig. (2) (from Hayashi's work). The periods spent in each 
region of the graph are marked on the curves. Such graphs 
serve as a basis to the study of nucleosynthesis of minor 
products in stars. 

We note that in a star of medium mass (4 M Q ) , hydrogen 
burning in the core and in the shell takes place at very 
similar values of p and T. In bigger stars the shell burns 
at a quite larger temperature and slightly larger density. 

Conditions in the contracting core are interesting for a 
study of the pre-helium-burning stage. 

The nucleosynthetic effect of the hydrogen burning stage 

is the transformation of H, He, Li, Be, B into He^ and the 

14 

transformation of C N 0 F (and maybe Ne) in mostly N (~95%) 


but also C 12 (~4%) and C 13 (~1%) . In fact at the end of 
H-burning these ratios are not very much temperature dependent 
(Caughlan and Fowler 1963) and are characterized by the 
equalities 


n 12 P 12,l n 13 P 13,l n l4 P 14,l n l5 P 15,l 


where P 12 ^ is the probability of the reaction C 12 (p,y)N 13 , 

12 

and n^ 2 is the number of C per unit volume. These equalities 

will be of importance to us later. In Pop II stars the 

. 14 -3 

fractional mass of N created this way is ~ 10 , while it 

is ~ 10~^ in Pop I stars . 


HELIUM-BURNING- REACTIONS 


C 13 Substaqe. C 13 (a,n)Q 16 . Q = 2.215 MeV . 

During the contraction phase the temperature becomes high 

13 16 

enough for the reaction C (a,n)0 to take place. The rate 
is approximately given by (Caughlan and Fowler, 1964) 


!°g (Pi 3 ,4 / PX 4 f 14 4 ) = 17.1 - 2/3 log Tg - 30.2 T 


-1/3 


8 


iog (e 13f4 /X 13 ) = log P^ + 17.2 


(13) 


(14) 

(15) 
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The rate is plotted in Fig. (3) (without electron screening, 
an important factor in small stars) . The concentration of 

varies from ~ 10” ^ (Pop II) to 10”^ (Pop I) . The 

13 

energy released in the entire fusion of C is 171 X^ KeV 
per nucleon in the gas. Even in a Pop. I star and even if 
the core material is highly degenerated (so that the nuclear 
energy release will heat the core instead of expanding it) 
the rise in temperature would be less than one million degrees. 
Consequently this reaction seems to have rather little in- 
fluence on the course of stellar evolution. 

Alan Liebman has used a sequence of models by Schwarzs- 
child (1962) to investigate the onset and importance of this 
reaction (and of others to be described shortly) in a Pop. II 
star of 1.3 M Q on its way to the helium burning phase. In 
Fig. (4) the energetic evolution of this star is described. 

L/M is the overall energy generation rate while L core / M core 
follows the contracting helium core itself (here most of the 
energy comes from hydrogen burning in a shell) . In the model 
C (and N X4 ) burning are neglected. From the graph it is 
clear that even the largest reasonable amount of (X = 10“^) 
cannot play a significant role. 


14 14 18 18 j. IP 

N Substage , N (g.y)F . Q = 4.404 MeV, F (B , v)0 ,Q = 1.677 MeV . 

14 

Later in the contraction phase N becomes combustible. 

New data on this rate are analyzed in Appendix A3. In the 
range 0.4 < T < 1.7 the rate is given by 

O 

l0 % P 14.4 / "Vl4.4 = * 3 ‘ 8 ' 12 ' 45 7 T 8 * 3/2 109 T 8 (16) 


1o 9 6 14, 4 7 X 14 ’ 17 ‘ 7 + 1o 9 P 14,4 ’ 

n = 1 n = [28.7/T - 1 . 5] (unscreened) 
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(17) 


As this reaction will turn out to be important in small 
stars (high densities and low temperatures) the electron 
screening factor (f^ ^) will be large. In many cases one will 
meet the conditions of strong screening (Salpeter 1954) . 

The energy released by the fusion of N is 434 X keV 

14 14 

per nucleon in the gas. In a Pop. I star (X^ ~ 10 ”^), with 
degenerate core, this could correspond to a rise of several 
tens of millions of degrees. In Fig. (4) the effect is plotted 
in the sequence of models of Dr. Schwarzschild (1962) (assuming 
naively that this extra source of energy does not react on the 
core) . Clearly this reaction could trigger the helium flash 
(the lifting of degeneracy) . In smaller stars it could even 
achieve the entire flash by itself. The flash would best be 
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called the nitrogen-helium flash. 

Hayashi has shown that stars with M < 0.53 (Pop. II) 
or M < 0.42 M 0 (Pop. I) cannot lift their temperatures high 
enough to allow the occurence of the 3He 4 -* C reaction, 
hence to qualify for the helium burning stage. The effect 
of the N^ 4 (He 4 , yjF 1 ® reaction (neglected by Hayashi) should 
lower these limits, especially in the case of Pop. I stars. 
This will increase the differential mass range of white dwarfs 
of both populations. 


Plasma Neutrinos 

Around T =1 the emission of neutrinos by the plasma in 
8 

stellar interior becomes important. Detailed calculations 
have been made by Inman and Ruderman (1964) . The results are 
shown in Fig. (5). 

One important parameter is the ratio (cp) of the plasma 
frequency u) o to the thermal energy 

-4 1/2 p . 2/3-i —1/4 

cp = TT oo q /kT = 3 ' - 3 ^ x - - — (P/n e ) I 1 + 1.0 x 10 ( P /ji e ) J 


(18) 


Here p is the mean number of nucleon per electron (p ~ 2) . 
e e 

The line cp = 1 is plotted in Fig. (1). In the region cp > 1 
the density and temperature exponents are given by 


m ~ 2.75 - cp/2 , n ~ 1.5 + cp 


(19) 
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The temperature exponent is always much smaller than the 
exponents associated with nuclear reactions. 

In Fig. (1) the plasma neutrinos are seen to be most 
intense around cp ~ 5.5 (summit of the ridge). In Fig. (4) the 
energy dissipation by plasma neutrinos in the core of a 
1.3 star during the pre-helium phase is plotted. 

During most of the contraction phase is about ten 
times smaller than the gravitational energy generation of the 
contracting core, but it reaches about one third of this value 
at the end of the contraction, after the nitrogen flash. It 
seems that the nitrogen flash comes just in time to prevent 
this process from influencing appreciably the course of the 
evolution. 

, 4 12 x 

Helium Burning Stage. (He -» C ) 

4 12 

Around T =1 the reaction 3He -» C becomes an important 
8 

source of energy. The rate can be computed by methods of 

statistical mechanics or by ordinary reaction rate theory. 

The second way seems to give a rate twice as small as the 

first one, if one overlooks the fact that the Breit-Wigner 

cross section formula for identical particles (in allowed 

states) is twice as large as for non-identical particles (the 

2 2 

geometrical factor is 2rrk instead of tt& ) . When this is taken 
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into account the results agree 


, < 3. ,3/2 6 -Q/kT 

( n i2 * /n 4 ) = 3 \ a e 


where * is the number of excited C nuclei in a state at an energy 

Q above the mass of 3Hei and X. = h/(2irM kT)^^. 

a a 


The properties of the second excited level in C have been 

reanalysed by the Brookhaven and the Cal. Tech, groups. The level is 

an 0 + with T = 2. 4 + 1. 5 meY, and Q = 372 t 4 keV . The reaction goes 

resonantly through this level in the range 0.80 < Tg < 6. 

12 

The rate of creation of C 


1°9 P 3a-C 12 /(PX^) 2 f 


4' 3a-C 


= -(6.72 + 0.25)-18.9 T" -3 log T Q . 


and the rate of energy generation is 


!°g c = f 11 - 52 ± 0 . 25 ]+ log[f 


p ^ 

3<x-C M 4 


1 -3 log T - 18.9 T" 
J 8 J 


The uncertainty on the rate comes partly from T 1 and from Q. 
The exponent m = 2 while n (unscreened) ia given by: 

/ 

n = 3 + 44.5/Tg, a very high value indeed. 


The energy generation rate during helium bu rning also depends 
on the subsequent reaction (& , y ) * 'This is still a weak 

point in our theory as its rate is poorly known. It can be written as 

108 [ P l 2 . 4 / P X 4 f 12, 4 ] = 10 ' 7 + 108 (S a> ~ 30 ' 05 T 8 " Z log T 8 (23 > 
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2 12 4 

The parameter ©^ represents the fraction of the C + He state in the 

7. 112 MeV level in O^. Nothing is known about its value except that 
2 

©^ < 1 and most likely > 0.01. It could in principle be determined by. 

experiments but so far it has not been. Some authors have discussed 

2 

its possible value in terms of alpha particle models (0 — 1) or collect - 

CL 

2 

ive models (@ 7® or possibly 0.024). At the present time it seems 

2 

safer to use an intermediate value such as % =0. 1. The idea here is not 

a 

to use an "average value" but to minimize the possible damage made on 
further theoretical development by an eventual experimental determination 
of this quantity. 


On the 0 16 (He 4 ,Y)Ne 20 reaction (Q = 4.730 MeV), some 
very recent data is discussed in Appendix A 4. The rate 
becomes (in the range 2 < Tg < 8) 

log P 16,4 7 PX 4 f 1 6,4 = < 2 - 1 1 °- 2 ) x 1q2 T 8 " 3/2 10" 45,9/T 8 (24) 

it is three time smaller than the rate given e.g. in (Reeves 1964). It 

12 4 

is so much smaller than the rate of C 4 He (for any reasonable value 

2 

of © ) that except for very heavy stars M > 30 the helium process 
never goes past the former reaction. Consequently we can write the 
total energy generation rate as: 


e = e 3H 4 +Ci2 
e 


y - Cl I 1 ° 12 - 4 

a V 3Q 3He 4 * < C 12 


771 q 16 > ) s C 1 + ° - 33 
X 4 




(25) 


(26) 
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l°g(q 16 p)=[lc<P 12i4 /p X 4 £ 12 ^-log 12 /(Px/f (27) 

ja — v-> 3 q— C 

12 

The term q 12 *- s the rate °* formation of C equ. (2 1). 

Here (q, p ) is independent of the density and of electron screen- 
1 6 

2 

ing effects. It contains the term 9^ and hence is very poorly known (about 

a factor of ten). This usually has rather little effect^ on the total energy 

generation rate since ¥ stays below two or three during most of the 

Helium burning phase. However Deinzer and Salpeter (1964) have shown 

that at the end of Helium burning (X^ < 0. 1) ¥ will become very large, 

2 

and there the uncertainty on •O plays a major role. 

To summarize: for stars with M < 30 the energy generation 

rate in Helium burning phase is given by eqs (25), (26), (27) and (22). 

The uncertahty in the rate is about a factor of two for £ 0. 3. For 

X ^0.1 the uncertainty could reach a value of order ten. This uncertain- 
a 

ty may in turn deeply effect the determination of the trajectory of the 
star in the fcL EL chstgxam 

Helium reactions in massive stars (He^ , y )Mg^ 

0 = 9.314 MeV); Mg 24 (He 4 , y ) Si 28 , (Q = 9.986 MeV);: Si 28 (He 4 , Y ) S 32 , 
0=6. 946 MeV) . 

In stars with M^- 30 M© these reactions start contributing 
significantly to the energy generation rates. 

Methods to evaluate these reactions and the uncertainties 
attached to them are described in Appendix A. We obtain 

log (P 20,4 / P X 4 *20,4* = 19,7 “ 43.75 T g “ 1/3 -0. 09 T ft 2/3 - 2/3 log T r (28) 


8 


8 


log(P 24,4 / P X 4 f 24,4 } = 20 * 7 “ 49JB9T“ 1/3 - 0. 10Tg 2/3 -2/3 log T g (29) 


20 


4 


log(P 28,4 / P X 4 ^8,^ = 22 '° " 3 logT 8 ” 55 * 67T 8 1/3 ~ °- 12T 8 2/3 (30) 


In this case Y a in equ. (26) becomes 

—2 

'\l/> - , + r X 4i X 4 Q 4(i+1) 

U3 ~ 1 °:.4 „12 


q 4(i+l) 


(31) 


3He - C 


where 


log [q 4(1+1) p] = 1<- P 4i,4/ P X 4 f 4i.4)- lo s(P 12 /(P £ ) 

I • 3a-> C 

Here , clearly, i = 4 represents O 1 ^*, etc. With the proper Q values ¥ 


(32) 


As given by 



a 4- 

0.162X 

X 

q + 0.256X X 

12 

x6 

4 

16 

20 4 20 


q + 

0.136X, 


q__ 

24 

h 28 

4 

28 

32 . 


Nucleosynthesis from He lium -Burning -Reactions 

The regions of the p -T plane where Helium burning reactions 
are likely to be found (either in the core or in a shell durii^j the carbon 
burning stage) are displayed Fig. ( 6) , again from Hayashi’s work. (It 
should be remembered that neutrino emission has there been neglected. 

Its presence would affect the conditions of" helium shell burning but 
not the conditions of helium core burning . ) 

To evaluate the nucleo synthetic effect of Helium the isotopic 
evolution has been followed till X q - 0.01 assuming fixed temperature 
and density Fig. ( 7 ) . The curves are iso -abundance curves of various 
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elements. For instance the words > 50% in a certain region mean 

that in the corresponding range of density and temperature the final pro- 
duct is mostly O 1 ^. We noticed that Ne^ is almost by-passed; it appears 
only in a small region where it reaches at best a value of 10%. (By 
using extreme values for the rates a value of 20% could be obtained). 

Salpeter and Deinzer (1964) and Hayashi (1962) have followed 
in detail, the isotopic evolution during Helium burning in the core for 
various masses. Salpeter and Deinzer start with initially homogeneous 
Helium stars, a simplification which may be of importance in assessing 

the accuracy of the predicted isotopic results. For the 0^(He^,\)Ne^ 

2 2 

rate they use successively 0 = 0. 1 , 0. 4, 1 . (Fig. ( 7 ) uses 9 =0.1.) 

a a 

i Z 1 6 

Their results are shown in Fig. (8 ). The final abundance of C and O 

2 

depends very much on 0^ but not the final abundance of heavier material. 
The isotopes Ne^ and Mg^ are essentially absent of stars with M<30M^. 
Even for larger masses Ne^ reaches at best 20 °/ 0 (these authors used 
older data which gave an 0^(He^,\)Ne^^ three times larger than the rate 
quoted here) . 

1 Z 1 6 

Hayashi 1 s results on C and O in the Helium burning core 

20 

are quite similar. He does not consider the generation of Ne . Hayashi 1 s 

12 -16 

group has also studied the formation of C and O in a helium burn- 
ing shell during carbon burning. For stars of 4 and 15. 6 the tempe- 
rature in the shell is very close to Tg = 2, the density ^800 /the fraction- 
al carbon production is somewhat less than during helium core burning 

20 

and the production of Ne (evaluated from their models) is at best 
a few percent. 

The purpose of discussing these elements was to determine 
potential fuels for the following stages of evolution. It is of interest 
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however to attempt a comparison between these results and 
observational evidence on cosmic abundances . For such a 
comparison to be meaningful one would need a detailed theory 
of the restitution mechanisms of the evolved material to 
interstellar gas. Preliminary investigations have already 
been made on this matter. (Salpeter 1959) (Schmidt 1959) . 

These authors make the assumption that stars with M > O.7M 0 
restitute an amount AM = (M-O.7M 0 ) to the interstellar gas at 
the end of their lives. Assuming further that the relative 
rate of star formation as a function of stellar mass has not 
changed since the beginning of the galaxy (the total rate has 
changed!) one finds that most of the gas came from stars with 
2 M g < M < 5 M . (Stars with M > 10 M Q contribute virtually 
nothing.) The outer envelope, somehow expelled from such stars 
will contain, amongst other things, the product of helium 
burning reactions (mostly, presumably, from helium burning 
shells during later stages of evolution) . The physical condi- 
tions in these shells still have to be computed, taking into 
account the neutrino emission processes. However from Figure (7) 
we can safely guess that the neon production will always be small. 

Observationally, using the abundance of C, 0, Mg in the sun 
(Goldberg, Mueller and Aller,1962) and the abundance of Ne^ 
from Suess-Urey "cosmic abundances" (Si is the standard) it 


appears that of the group C,0,Ne,Mg, Si, about 20% belongs to 
carbon, 50% to oxygen, 20% to neon, and a few percent to 
magnesium and silicon. In view of the previous discussion the 
abundance of C and O are understandable but the abundance of 
Ne is very high. Since neon is also a product of carbon burning, 
most of it may come from there. Magnesium and silicon most 
likely come from the carbon burning stage, or from more ad- 
vanced stages of stellar evolution. 

Solar cosmic rays have recently been analysed by Biswas 
and Fichtel (1963) . The analysis gives strong support to the 
view that, except for proton, the accelerating mechanisms do not 
disturb the original relative abundances. They find indeed 
the ratio of C/O in cosmic rays is very similar to the solar 
atmosphere C/0 ratio. They have measured Ne^® (not spectro- 
scopically detectable in the sun) and found it to be smaller 
than the above quoted value by a factor of about three. With 
their value one obtains: X c *= 0.26 X Q = 0.58 X Ne = 0.09 


X Mg ' °- 04 S si * °- 04 - 


The cosmic "normal" abundance of Ne^® has been evaluated 
from analyses of B stars, which according to Dr. Stromgren are 
very uncertain. The planetary nebulae typically have ratios 
of neon to oxygen which are about 10 times smaller than "normal" 


(Aller, private communication) . It does not appear unreasonable 



to argue that if there are any settings in our universe where 
we may hope to see the products of helium burning, unadulterated 
by nuclear processes associated with further stellar evolution, 
it is indeed in the gases which form the bulk of a planetary 
nebula . 

Heavy elements build up in the helium burning phase . 

The presence in nature of elements with A > 56 is 

attributed (BBFH 1957) to the turning on at a certain period 

(or periods) of stellar evolution of a rather large flux of 

neutrons. In particular, the existence of isotopes in the 

bottom of the valley of nuclear stability is proof of the action 

of neutron flux with characteristics such that the neutron 

capture time was much longer than beta decay lifetime of the 

unstable isotopes lining on both sides the bottom of the valley, 

and such that several tens of neutrons were absorbed per each 

seed nuclei (usually the metals, Co-Fe,Ni*). This process is 

called the s process, and the elements thus produced are called 

13 1 fi 

the s elements. In this respect the reaction C (a,n) 0 
may be of importance since it generates neutrons, and hence may 
be a source of heavy element build-up. The thus burned may 

have been left from H burning (as discussed before, the frac- 
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*For more details see Clayton et al(1961), Seeger, Fowler and 
Clayton (to be published) . 



tional mass of the CNO group to be found in the form of C 

is about 1%) or it may have been formed by mixing some protons 

to the helium burning core (Greenstein and Wallerstein, 1964) . 

13 . 

The effect of the originally present C is negligible for two 

reasons. First, even if every neutron were to be used for 

heavy element build-up, at best one neutron would be absorbed 

per metal nucleus. Second, there is the question of the 

poison, about which however one must be careful. From Figure 

13 

(3) it is evident that the atoms of C will always burn before 
14 

the atoms of N , hence in the presence of these atoms. It is 

correct to say that through the N^(n,p) reaction and its 

14 

large cross section N will act as a neutron poison. But if 

12 

the proton thus produced is then captured by a C atom, we 

shall have regeneration of the lost neutron by the reactions 

C 12 (p, y) N'*' 3 - C 13 (He^,n) 0* 6 . In other words C 12 will act 

14 

as an antidote to N poison, if it is abundant enough to 

14 

capture a large fraction of the protons from the N (n,p) 

reaction. Will that be our case? Clearly no, since in view 

of the hydrogen-burning-equilibrium-abundance ratio prevailing 

in the core (see eqn. 13), the proton has at best a 25% chance 

12 

of being captured by a C neutron. Hence the formation of 

13 

heavy elements from the original C is to be ruled out. 


26 


After the onset of helium burning the abundance of C x * 

will increase sharply and the situation would be reversed 

if C 13 and had managed to survive (which will never be 

the case) or if they had been produced by admixture of some 

12 

hydrogen from the envelope. In this case C is a good 
14 

antidote to N poisoning and heavy element build up is a 
possibility. Such mixing has been considered recently by 
Wallerstein and Greenstein (1964) in an attempt to explain 
the peculiar abundances found in CH stars. 

In these stars carbon and the product of neutron capture 
process on a slow time scale (Ba, La, Ce, Nd) are enhanced 
considerably while oxygen is normal. The enhancement of C 

may imply that the neutron capture processes took place after 

4 X2 16 

the onset of the 3He -* C reaction. The absence of 0 

enhancement may imply that these processes took place rather 

4 12 

shortly after the onset of the 3He -* C , and most likely at 
rather low temperature. As we shall discuss later, this 
gives strong support to the mixing mechanism advocated by 
Greenstein and Wallerstein. 
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Neutrons can also be generated by the sequence 

0 18 (a, n)Ne 21 Q = -700 keV 


N^ 4 ( a , y) F^ 8 (e+ v) 0-*- 8 




O^- 8 (a, y)Ne^2 , 


Ne^ (ajnjMg^ 5 Q = -482 keV 


,25 


[ Ne^^ (a, v)Mg^ 8 


Because if its high threshold energy the branching 

toward Ne^l is very difficult to achieve. It would occur if 
14 

N atoms were to be suddenly brought to temperatures Tg ^ 4 

even for only a matter of seconds. Such conditions may be 

met during the helium flash. In Schwarz schi Id 1 s star (M = 1 3M^) 

T = 3.4 is attained for a second or so. Evidently more detailed 
8 

models will be needed before we can evaluate the 
importance of this mechanism in stars of smaller masses. 

For heavier stars the lower branch is followed. 

Two questions here: how far will the chain go before the 

supply of He 4 is exhausted; if the chain is carried to the 
end will the final products be Mg^6+ y or Mg^5+ n f To answer 
these questions we must remember that at the end of helium 
burning the main consumer of He 4 is . From the graph it 
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# 14 18 12 

^PP^srs that N uiid O wx 11 always burn be fors C f but not 

necessarily Ne 22 . If the process takes place at Tg < 2 

22 oo 

Ne will be left intact; at higher temperature Ne ^ bums 

2 5 

and from the graph again we get mostly Mg + n, with a small 

fraction of Mg 2 ®. With the abundances of discussed before, 

this chain of reaction can yield as much as 100 neutrons per 

seed nuclei, probably enough to meet the requirement of the 

s process. At still higher temperature (T >3) Mg25 itself 

8 

25 28 

is reacting by Mg (cL, n) Si to essentially double the output 

of neutrons. As a result we do expect sufficient neutron 

14 

generation from N isotopes if the helium burning stage 
reaches a temperature of at least T_ = 2 for a time sufficiently 

O 

long to bring about a fairly complete exhaustion of Ne 22 isotopes. 
Such temperatures are reached during the nitrogen-helium flash 
but for times which seem to be far too short (by factors of 
thousands, for instance in Schwarzschild 's models). Studies of the 
flash for different masses may give a different answer but 
at present this seems to be extremely unlikely. 

The proper combination of time and temperature are reached 
at the end of helium burning stages at least for rather massive 
stars (e.g., M = 7M^ as studied by Hoffmeister et al, 1964) . There 
however, we should expect the presence of s process elements 
such as Ba La Ce Nd to be accompanied by enhancement of C and 
of O. 
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To summarize: neutron will not be produced prior to the onset 

* 

14 

of helium burning. Small stars may get N generated neutrons 

TO 14 

during the helium flash; bigger stars do not get C or N gener- 
ated neutrons in the early days of helium burning. Neutrons may 
appear through admixture of hydrogen, producing (if the enriched 
material eventually reaches the stellar surface) enhancement of 
carbon, and of s elements, but not necessarily of oxygen. Finally 
generated neutrons will be produced at the very end of helium 

burning if the mass is large enough to bring the central temperature 

4 

above T 0 = 2 before the He is exhausted. From existing models, at 
8 

least a few solar masses are needed for this to occur. 

Photoneutrino Contraction Phase 

This name is probably appropriate since in the contraction 
phase between helium burning and carbon burning this process be- 
comes active and in fact, dissipates more energy than the photon 
emission processes. The switchover usually takes place in the 
range Tg = 3-5. 

The rate of energy generation is given by (Chiu, 1961) 

e v = 10 8 T g 8 /|a e erg/gm/sec (33) 

The term n e is the mean number of nucleons per electrons in 
the gas. Usually \ jl q - 2 (as for He^, C^, etc.). However 

at higher temperature electron positron pairs will be produced, 
which must be taken into account in evaluating u . 

The density and temperature exponents are m = 0, n = 8. 

Again this is a rather small temperature exponent, as compared 
with the helium burning (n ~ 30) and the following carbon 


burning (n ~ 30 also) . This contraction phase will end when 

the reaction C * + C will become fast enough to stop the 

contraction. Shall we have a carbon flash? According to 

Hayashi, Hoshi and Sugimoto (1962) , no carbon flash will take 

place if the photoneutrino process occurs. The core will 

quietly warm up till the carbon burning reaction equals the 

neutrino energy dissipation. 

Carbon Burning Reactions 
The Chalk River Group has recently studied very 

thoroughly (experimentally and theoretically) the reactions 

between C^ 2 and a few light nuclei (C^ 2 , N^ 4 , 0^ 6 ) . The 

elastic scattering data and the capture cross section data 

have been improved to about 20% accuracy. Optical model 

parameters have been obtained which match the data all through 

the range of energy considered. In the Appendix A8 , a method is 

described which allows a determination of the low energy part 

of the capture cross section and hence a determination of the 

thermonuclear reaction rate , most likely accurate to better 

than a factor of three. The probability of a reaction is 

given by 

m 

m 1/3 

log P/pX c = (26.4+0.5) - 36 + °r° . 70 =LL_ -2/3 log T Q ( 34 ) 

10 c ~ Tg y 

and the energy generation rate by 


log e c = 17.7 + log( p c+c X c ) erg/gm/sec 


(35) 



The range of temperature and density at which a star will 
stop its contraction to burn its carbon can be investigated 
by means of the following crude model. We assume first that 
the energy is almost entirely dissipated by neutrino emission. 
Several detailed models have shown this assumption to be valid. 
Then we find in the p.T plane the locus of all the points 
representing conditions at the center of the star where carbon 
burning energy generation (assuming pure carbon) equals the 
neutrino energy dissipation. For temperature neighboring 10 9 
we should include pair annihilation neutrinos (see Figure 1 ) 
which we shall discuss presently. This curve appears in 
Figure (10) under the label "first model". Correction has then 
been made to allow for the fact that neutrino emission takes 
place in a bigger volume than nuclear energy generation. This 
is the "corrected model" curve. 

Hayashi has published models of stars in the carbon 
burning phase for which he neglects however the neutrino energy 
dissipation. He considers stars of 4M Q and 15.6M Q , and sets 
at O.7M 0 the lower mass limit for a star to reach the carbon 
burning stage. In Figure (6 ) the evolution of the central 
density and temperature (together with the helium burning shells) 
are pictured. The H.R. diagram of the cluster h and x Persei 
can be used to test the validity of the model. 
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* There , two well defined types of giant stars appear; a group 
called early-type supergiants 22 stars) and a group called 
late type supergiants ( ~ 15 stars) . In between lies the 
Hertzsprung gap. Hayashi's evolutionary track after leaving 
the main sequence goes slowly through the first group, rapidly 
through the Hertzsprung gap and then slowly again through the 
second group. We expect the ratio of the period spent in each 
group to be similar to the ratio of the number of stars in 
each group. The ratio of the periods is about one in Hayashi's 
model (neglecting neutrino emission) . Hayashi has calculated 
(although crudely) the reduction in the carbon stage period 
resulting from neutrino emission. He finds that the number 
of stars in the corresponding group should be down by a factor 
of about ten. Detailed models of carbon burning stars with 
neutrino emission would be desirable. If similar results are 
found, such evidence would work against the assumption of 
direct coupling between electrons and neutrinos (although the 
statistics here may not be really relevant) . 

Nucleosynthesis during carbon burning bhase . 

The reaction C^ 2 + C^- 2 produces Ne 2 ® + He 4 and Na 22 + H 

in roughly equal amounts. There is also a weak endothermic 

23 

branching toward Mg n which we shall discuss later. 
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The abundance equations of carbon burning have been 

integrated for various carbon and oxygen initial abundances 

(X^ 2 ^ + x i6^ = 1) for different sets of densities and 

temperature by Cameron (1959) ; Reeves and Salpeter (1959) 

and Tsuda (1963) . The results are rather insensitive to 

the choice of p and T and to the (slightly different) 
c c 

choice of rate made by each author. They can be qualitatively 
described in the following way: (i denotes initial abundances 
and f denotes final abundances) (a) X^ 6 remains the same 
(X lS ^- ~ X 16 ^) ; (b) x 2Q f is always less than 0.35; for X^ 1 <0.40 

X 2 0 1 /X 12 1 p 0.6, -for X^ 1 >0.40 : X 2Q f -0.30; (c) X^/X^ ^ 0 .40 

If the carbon reaction takes place at relative energy 

. 23 

larger than 2 . 6 MeV the branching Mg +n occurs about 5% of the 
time. (Bromley et al 1960) When the stellar central tem- 

perature is less than T = 0.8 the neutron flux thus generated 
is negligible (see Appendix B) . 

On the other hand at Tg<0.75 substantial fluxes of neutrons 
become available through the following set of reactions 
initiated by the protons and He^ emitted in the main branching : 

C^ (p» y)N^ (e + , v)C^ (a,n)0^. Calculations have shown that the 
12 13 

C (p,y)N reactions occur mostly in the early days of the 
stage, and consequently in the lower temperature range of carbon 
burning. This way neutron fluxes capable of inducing about 50 
neutron capture per metal nuclei could be released. 


This source of neutrons is quenched at T e >7.5 as then the 

o 

13 . 12 

isotope N photodxsxntegrates into C +p before it has time 

13 

to beta decay to C 

The introduction of pair annihilation neutrino processes 
has recently called for a reevaluation of the temperature at 
which the carbon burning stage will take place. 

A combination of Hayashi's model (inhomogeneous models 
without neutrino emission) and Salpeter and Deinzer models 
(homogeneous models with neutrino emission) makes it clear 
that the carbon burning temperature ranges from Tg = 0.8 to 
about 1.1. In the table the number of neutrons emitted per 
metal atom, assuming - 1/2 X m = 5x10“^ is given 

T g 0.6 0.7 0.8 0.85 0.9 1.0 1.1 1.2 

No. of 

neutrons/metal 0.1 0.6 4 9 13 20 40 50 

12 12 

Clearly the C + C burning stage provides an excellent source 
of neutrons and s element processing. Here the earmark would 
be probably the simultaneous enhancement of s elements and of 

23 a j 

Na (especially when compared to Mg ) . This statement cannot 
be proven simply. Qualitative analyses have shown it to be 
true in most cases of interest. In a naive way, however, one 
sees that each neutron absorbed leaves behind an Mg 2 ^ which 
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23 

quickly decays to Na therefore increasing the concentration 

of this isotope. In real life one must still consider further 

23 

reaction between all constituents of the gas. Na usually 
turns out to be still enhanced (as compared to the case where 
the branching ratio is shut) but not nearly as much . Detailed 
calculations would be needed to settle this matter. 

From the previous paragraphs, s processing taking place 
during the carbon burning stage should also lead to enhancement 
of C (if the material is extracted before the end of the 
stage) 0^ 6 , Ne^ and Mg^ together with enhancement of s 
elements . 

Pair-annihilation-neutrino contraction phase . 

In the legion T =6 the dominant emitting neutrino process 
8 

passes from photoneutrinos to pair annihilation neutrinos 
(see Figs. (1) and (9)). For a nondegenerate nonrelativistic 
gas the energy generation rate is given by 

log e pair = 18.7 - 5.15/T g + 3 log T g - log p 

with exponents n = 11.9/T +3, m = -1. 

9 

Neon-photodisinteqration (flash ?) 

At Tg ~ 1.25 the contraction will be perturbed by the 

20 

photodisintegration of Neon . The lifetime of an atom of 
20 

Ne against photodisintegration is given by log^ 0 t = 


(36) 


36 _ 


, -(12.15 + 0.05) + 28.4/T . Computation (Tsuda 1963) has shown 

that one half of the destroyed Neon becomes 0^ and the other 

24 20 24 

half becomes Mg (through Ne zu (a,y) Mg 4 ) . The energy yield 
is about 2.3 MeV per initial neon atom. Hence 

log^ e = 29.2 - 28.4/T^ erg/gm/sec. ^ 


As discussed in the previous section the core at this 

20 

coment contains about 25% Ne (in fractional mass) . The photo- 
disintegration will yield ~ 600 KeV for each nucleus in the gas 
or ~ 30 KeV for each particle in the gas. In a partially 
degenerated core we may have a neon flash bringing us all the 
way up to the oxygen burning stage similar to the flash. 

For large stellar masses the Ne^° photodisintegration will at 
best slow down the contraction towards the next stage. 

At the end of this period, the isotopic abundance is 

^ • 

roughly as follows: (X^ 2 + X 16 = ^ describes the results of 

the previous helium stage). (a) ~ X * + 0.15; 

(b) X 20 * 0 ; (c) X 24 * 0.4 X 12 + 0.15; (d) X 28 = 2 to 4% . 


Oxygen Burning Reactions 

As the core temperature approaches Tg = 1.3 ~ 1.4 reaction 
between oxygen nuclei start releasing energy. Here experimental 
data is rather lacking. Only the O^ (O^®, a)Mg^ has been 
measured (Almqvist 1960) and the data extends only slightly 
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below the Coulomb barrier. More work is being done on that re- 
action but is not available as of now To an accuracy no better 
than a factor of ten the following equation can be extracted from 
the data 

log P /PX Q = 40.5 - 59.02(1+0.14 T g ) 1/3 -2/3 log T g (38) 

with 

log e Q = 17.8 + log (P X Q ) . 

In Figure (12) the temperature and densities at which oxygen 
and neon are expected to burn are described by using the crude 
model described in the section on carbon burning. The nucleo- 
synthetic effect of the oxygen burning stage has been studied 

by Tsuda (1963) and Cameron (1959) . The main outcome are the 

24 32 24 28 32 

elements from Mg to S . The elements Mg , Si , and S 

are not really preponderant over in-between isotopes (as is the 

case in cosmic abundances) . 

After the fusion of oxygen the emission of neutrinos becomes 
so strong that nuclear energy generation most likely never 
succeeds in halting the contraction. Hence we do not expect 
any other burning stage. Nuclear reactions leading to the 
equilibrium process have been considered by Hoyle and Fowler 
(1964) and Tsuda and Tsuji (1963) . They will not be discussed 
here . 
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APPENDIX A 


Thermonuclear reaction rates in the range T = 1 to 50 

For some (a,y) (a,n) (p, y) reactions, e.g.. Mg 24 (a, y) Si 28 

the level scheme of the compound nucleus has been studied 
experimentally in a rather wide range of energy and the widths 
^a'^p'^y and t ^ ie statistical factor uu = (2J+1)/ (2I t +l) (21^+1) 
have been determined for all the resonances E r in that range. 

In such a case a lower limit to the rate can be obtained by 
adding the effects of each resonance 

P/PX = E Pi/pXi 
all 

resonances 

r r n 

log Pi/pxx = 12.69 - 3/2 log p - log A. + log Q* -? n ) ( A l) 

- 3/2 log T g - SO^E^Tg . 

Here pi is the reduced mass, and r c , T n the widths for the 
charged particle and neutral particle reactions, respectively. 

This rate is .a lower limit because of the contribution of reso- 
nances lying outside of the experimentally explored range. 

Because of experimental difficulties the (a,y) reactions have not 
been brought below 1 MeV and the resonances in that range are 
expected to dominate the rate at T < 3 or so. In Fig. (15) the 

O 

• oa 

position of the lowest measured resonances in Si are plotted 

against the background of the Gamow peak to illustrate the 

situation. 
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The total capture rate can also be evaluated by a combi- 
nation of the information on individual resonances with data 
obtained from optical models of the nuclei. Here we take 
advantage of the experimentally discovered regularity of the 
optical model parameters, regularity best expressed in terms 
of the strength functions. In opposition to the case of neutron 
reactions, alpha reactions have so far failed to show structure 
in the individual partial wave (l) strength functions. The sum 
of these strength functions appears to be energy independent, 
at least in ranges wide enough to be of use in the evaluation 
of thermonuclear reaction rates. 

Through the optical model a computation of the capture 
cross section can be made. For alpha particle reactions, the 
nuclear potential is chosen to be of the form 

V (r) = - (V Q + iW c )/ {1 + exp [(r-R 0 )/ a ]} . (A2) 

The best choice of parameters appears to be 

V 0 = 50 MeV Wo = +10 MeV 

Ro =(1.25 A t 1//3 + 1.6) f. a = 0.50 f. 

(A is the mass number of the target nucleus*) (Vogt, private 

communication) 

For proton reactions the real potential has the same form, 

1/3 

with V q = +55 MeV, R q = 1.25A t , a = 0.5f/ The imaginary po- 

tential has a gaussian shape centered on the nuclear surface with 
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W = +4 MeV. Calculations have been made by Chalk River workers 

y 

on a Bendix G-20. 

The results of these computations are then fitted by a formula 
of the form: (We use o in barns, and E in MeV) 

o = S e -2nn e -gE _ S -0.98935 (Z 1 Z 2 y 1/2 )/E 1 / 2 e -gE (A3) 


where g is left as a parameter determined in such a way as to 
make S energy independent up to the highest possible energy. 
Usually one can find a g such that S varies by less than 10% 
when the energy is varied from zero to about 60% of the Coulomb 
barrier energy (B = z 1 z 2 e / R ) • Then S becomes a measure of the 
sum of the strength functions. If the computation had been made 
with the choice of a square well "black nucleus" the factor S 
would have taken the form (see Feshbach 1953) 


S s.w= 2 -p 2 [?^<2M)G> 2 ]^ 


(A4> 


2 2 

The first factor can be written down as 2tt X /E = 4.126/V, |j being 

the reduced mass. The term in the brackets represents the 

penetrabilities for all t waves for a square nucleus of radius R. 

2 

is the irregular Coulomb wave function. § , an almost energy 

independent parameter, is defined as 

§ 0 2 S Y e2TTT| /TTG 0 2 = 2y 1/2 e 4yl/2 , y = 2R b /R; R b = [hVuZ^e 2 ] . 


The last term in eqn. (A3) is the value of the strength function 
for a black nucleus; K ~ 1.4f. is the wave number inside the 
nucleus. 
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At energies much smaller than the Coulomb barrier and 


for y 


!/2 > 


4 (as will be the case for our work on alpha capture 


12 12 22 
reaction and C + C ) the factor £(2-t,+l)g Q /g becomes 

independent of energy. Its value is plotted in Fig. (13). 

Numerically the value of S is closely given by 

S / w 

!°g S s w = 0.1 + 1/2 + 0 • 458 [ 55^2 (iR] 1 ^ 2 


(A5) 


The results of optical model calculations for alpha capture by 

12 32 

all the stable isotopes from C to S have been used to ob- 
tain the values of S Q m (optical model). All the values were 

very well represented by putting in the square well formula eqn. 

1/3 

(A5) the choice R = [1.50 A. + 2.0]f. (Fig. 14) 

s 1 

In the square well formalism, the value of g is 


3 

3/2 2 rt^R -i 1/2 _i 

a = y 3/2 /6n = 0.122 MeV x 

9 s,w Lz^J 


(A6) 


with R in f, i a in a.m.u. 

The values of g obtained from the optical model analysis 
have been plotted in Fig. (14) and tabulated in Table 1. It 
will be noted that they decrease somewhat rapidly with A^ the 
mass of the target. A sort of effective radius can be obtained 
by inverting the formula 


= 4.07(fi !2)' 

V n y 


Z^l/3 2/3 


f . 


(A7) 


This time the radius R needed to reproduce g was of the form 


R g = 0.7 A x 1/3 + 4.0. 
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In a forthcoming paper (Vogt and Reeves, to be published) 
an analysis of the meaning of these parameters in terms of the 
structure of the nuclear well will be made. 

Once S and g are found, the thermonuclear reaction rate 
becomes (S in MeV barns) 


1/3 1/3 

log P/P\ i = 10.56 + log Q-l -) -2/3 log Tg + log S 


- 3.975 ^ (1 + 8.62 x 10” 3 

T 8 


g Tg) 


1/3 


(A8) 


In Table (1) the values of S, g and also the numerical 
constants necessary to compute the thermonuclear reaction rate 
are given (some of them have been slightly modified by the 
results of other techniques). To identify these constants we 
rewrite the eqn. (A8) in the form 


log P/px = M - 2/3 log T 0 - N T 0 

4 o o 


- PT 


8 


(A9) 


Although the optical model usually gives a fairly good 
representation of the data, individual nuclei may at times 
depart from their expected behaviour. Then both S and g could 
influence the rate to any great extent. 

As mentioned before, from the experimental determination 
of the individual resonances an evaluation of S can be made 
directly. If the resonances do not overlap to the point where 
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log S (Mev b) 

g(Mev 1) 

M 

N 

P 

c 12 + He 4 

6.45 

.82 

16 . 6 

30.11 

.07 

C 13 + " 

6.55 

.83 

16.7 

30.25 

.07 

N 14 + " 

7.22 

.78 

17.4 

33.72 

.08 

N 15 + M 

7.29 

.79 

17.5 

33.88 

.08 

0 16 + " 

7.95 

.75 

18.1 

37.20 

.08 

o 17 + •• 

8.04 

.77 

18.2 

37.35 

.08 

o 18 + 

8.14 

.78 

18.4 

37.47 

.08 

p 19 + „ 

8.73 

.74 

18.9 

40.67 

.09 

Ne 20 + " 

9.44 

.72 

19. 7 

43.75 

.09 

Ne 21 + " 

9.39 

.72 

19.6 

43.88 

.09 

Ne 22 + " 

9.47 

.73 

19.7 

43.99 

.09 

Na 23 + " 

10.04 

.70 

20.3 

46.97 

.10 

Mg 24 + " 

10. 50 

.68 

20.7 

49. 89 

.10 

Mg 25 + " 

10.65 

.69 

20.9 

49. 97 

.10 

26 

Mg + " 

10.72 

.69 

21.0 

50.06 

.10 

27 

A1 + " 

11.25 

.67 

21.5 

52,91 

.10 

• 28 

Si + " 

11.^5 

.66 

22.0 

55.67 

.11 

Si 29 + " 

11.82 

.66 

22.1 

55. 75 

.11 

si 30 + " 

11.89 

.67 

22.1 

55.82 

.11 

31 

P + 11 

12.38 

.66 

22.7 

58. 56 

.11 

32 

S + " 

12.86 

.63 

23.1 

61. 16 

.11 


Table 1: Parameters for the computation of thermonuclear 

reaction rates. (The temperature unit is Tg = 10®° K) 
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interference phenomena become important, the cross section for 
individual resonances can be written in the Breit Wigner fashion 


a = tt* 2 « r c r n / (e -E r ) 2 + r 2 /4 r = r c + r n 


(A10) 


Averaging over resonances in an energy range AE 


<cr> = J adE / AE = E 2TT 2 t 2 T T ID /AEF = 
* c n 


(All) 


We define 


T = v ? 2 e “ gE e"’ 2TTT ' 
c,l y ct e 


2ttT| . 2 

= ye /iTG^ 


(A12) 


where y .is the reduced width and % an almost energy 
C / 

independent factor representing angular momentum effects (t) . 

If r > r all through the range AE, we have 
c n 


<Q> as S e -2Trn— gE 
E 


(A13) 


with 


»- ■ 


2ttVeE < ” . S 4 Yc xl > = [2tt 2 ^E] ? 0 2 (Y C/D ) = 


(A14) 


4.126 _ 2 

?o ^c/D^ 


Here D is the average distance oetween resonances of same (t) 

‘‘L / 

with orbital angular momentum (t) . The last term is a definition 
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2 

of % (y_,/D) based on the fact that the individual strength 

function < y^/D^ > are known to be practically energy 
independent. Using the experimental values of r a ( ex P) we 
calculate 


2 r (exp) 2ttti gE 

?o (Y c/d ) = S -L&- e Le y 

all 

res. in AE 


AE 


(A15) 


On the other hand, when T > T 

c n 

<ct> = 2ir 2 * 2 S — = 4 — 1 ? l 6 /Dl\l = £ (A16) 

AE AE |_i \ D / E E 


with C in MeV barns. The total rate involves the integral of 

<a>v over all energies. However the range of energies where 

T < F bring a small contribution to the rate if the Gamow 
c n 

energy E Q is above the energy (E*) for which F n = T . Then we 
write (if is Avogadro's number) 


P/PXi 




(2^kT) 


1/2 


Ce” E */kT 


(A17) 


or 


log P/p\ 1 = 10.00 - log - Vl log \s - 1/2 log T g + log C -50.4E*/T g (a18) 


The value of E* and C should really be obtained by trying to 
fit the eqn. (A18) with the values obtained in eqn. (Al) by 
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summing over all resonances. If the measurements of individual 
resonances are not available then we may use the value of S 
obtained from the optical model calculation to get an average 

r a 

value of and by equating it to r^/D obtain E*. That is 

S = 2tt 2 * 2 E § q 2 (y/D). 

TP = § Q 2 (Y/D)e~ 2TTT1 e*' gE = wr .for E = E* 

u V/D > 


(A19) 


We discuss next the validity of the optical model approach 
in the (low) temperature range where the widths of the Gamow peak 
become comparable to the average spacing between resonances. 

Except for a few nuclei, (p,y), or (a,y) reaction threshold are 
to be found in energy regions of the compound nucleus where the 
levels are fairly crowded. Typically the spacing of levels will 
be of the order of 100 keV but larger distances between two 
given levels will not be uncommonly found. 

The use of optical models for thermonuclear reaction rates 
clearly implies that at least a few resonances (or at least one 
resonance with a more or less average reduced width) are to be 
found inside the Gamow peak for any temperature under consideration. 
Otherwise the method may overestimate the rate by several orders 
of magnitudes . 
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The position and width of the Gamow peak are given by 
E q = 26.3 (Z^ Z 2 2 p) 1/3 Tg 2/3 keV 

AE q = 35 (Z x 2 Z 2 2 b) 1/6 ig 5/6 keV (A20) 

T' = [T 0 /l + 8.6 x 10" 3 g T Q ] g in MeV -1 

8 o a ) 

Here AE Q is the full width at half maximum, g is defined in A3. 

In the table the values of E Q and AE Q (in keV) are given 
for typical protons and alpha reactions. 

Z , Z 1,6 1,10 1,16 2,6 2,10 2,16 

1 2 

AE~ 64 75 90 100 120 140 

T-^2/3 
8 

E o 85 120 165 215 305 415 

m 5/ 6 
8 


Clearly in the region T g < 5 important departure from the 
assumption quoted in the previous paragraph may be expected. 

If the corresponding level scheme is experimentally unknown 
the determination of the rate becomes most uncertain. The rate 
obtained from optical model calculations becomes essentially 
an upper limit. A lower limit, an expected value and an estima- 
tion of the uncertainty can be be obtained by throwing in more 
nuclear statistics. 
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Experimentally the probability that two neighboring nuclear 
levels will be found at a distance D larger than a preassigned 
value D q is well represented by the Wigner distribution formula 

P ( D > D q ) = exp l-T/4 (D 0 /d ) 2 ] (A21) 


where D is the expectation value of D. The probability that two 
given resonances will be at D > 3D is about 0.001. For our 
purpose we restate this in the following way: our chances of 
finding one level at an energy E r within AE = Er - E 0 = 1.5 5 
of the Gamow peak E Q is 99.9%. We define as our upper limit to 
the rate of the value obtained if the level is at the center of 
the peak (the rate is then closely equal to our optical model rate) 
and our lower limit, the value is obtained if the level is at a 
distance E r = E Q + 1.50 D. The rate itself will be the geo- 
metrical means between these two values. The choice of 99.9% 
is arbitrary to an extent. It is however in the spirit of the 
determination of uncertainties by experimental physicists. 

Defining u = AEq/Eq, v = (3.0D/AE Q ) the ratio (f) of the lower 
limit to the upper limit is given by 


f = P (E = E n + 1.50D) / P(E = E ) 
x ° r o 

= exp [ (- fa*, (1 + uv + 7r 4 7I/2 -3) ] 


(A22) 
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Consequently, we have P max = P (optical model) as in eqn. (9) ; 

Pm-ir, = [P f] and we choose for the expected value and the 
mm u max 

uncertainty [log P + A log P] = [log P (opt. mod.) + 1/2 log f] 

+ 1/2 log f]. 

The reduced widths themselves have a distribution given by 

2 

P(V < Vo) - J ° -2 dx x = y/ 7 , 23 

(so that y = y + 0.9y with a probability of ~ 90%). The correct 
way of doing this problem would be to analyze the probability 
function of the ratio (y/D) . However this is somewhat compli- 
cated for our purpose. With our previous choice of D Q = 3D the 
chance of having two levels in the peak is itself 80% so that the 
fluctuation in the strength function averaged over the Gamow peak 
should probably not reach a value of more than three. A better 
estimate of log P should then be 

log P + A log P = ["[log P(opt.mod.) + 1/2 log f] + (0. 5+1/2 log f)J 

In cases considered here, where the optical model has any 
interest (O 18 (a, y) Ne 22 , Ne 20 (a, y)Mg 24 , Mg 24 (a, y) Si 28 etc.) the 
average level distance is about 100 KeV. Then 1/2 log f is 
about 0.8 at Tq = 1, 0.3 at T g = 2, 0.15 at T g = 3 and 0.1 at 
T 0 = 5. Since we shall rarely need these reactions at tempera- 

O 

ture T 0 < 3 , we have not considered this factor in the rate 

O 

given in Table 1. 


1 


(A24) 
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12 1 ft 

A-l C (q,y)0 The rate is discussed in detail in Cartledge, 

Thibaudeau and Reeves (1963) and also in Reeves (1964). 

13 16 

A- 2 C (g.n)O An estimate using an experimental measurement 

of the low energy non resonant contribution has been made by 
Caughlan and Fowler (1964). They obtain a value M = 17.1 while 
the optical model calculation yields 16.7. 


The properties of the excited levels of F J 


the vicinity of 4,4 MeV can be inferred through an analysis of 
the Ne 20 (d,a) and 0^(He^,p) reactions. (Enge and Wojtasek, 1959) 


In particular, a level with T = 1 should be identified by 

18 

the fact that the a + F outcome would be forbidden. 

The relative intensity of alpha groups from the Ne 20 + d 
reaction is given as 



Energy 

Re lative 


of the level 

Intensity 


MeV 



4.108 

125 

i 

4.218 

138 

! * 

4.350 

147 


4.400 

68 


4.651 

30 


4.741 

<7 


4.844 

70 
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Clearly the 4.741 level cannot participate in the reaction. It 

is most likely at T = 1 level. The 4.651 MeV level however is 

definitely active although its contribution is about four times 

smaller than the contributions of the other levels. Since the 

outcoming a have more than three MeV, the Coulomb effect is not 

important. The reduction probably reflects an unusually small 
2 

value of 0 , or a high angular momentum of the excited state, 

a 

or both. 

The reaction rate can be written as 

log P 14/4 /pX 4 = H.3 - 12.45/T 8 + log u)r a - 3/2 log T Q 

Following Brown (1962) we choose R = 5.6f for the radius of 

interaction. We obtain for uur of this level 

a 

log uuF^ = - 13.4 + log (w0 a ^) , 

2 

The average value of ) for levels of slightly higher 

energy is ~0.06. Here, because of the low experimental yield 
2 

we choose u)0 =0.02 and obtain 

a 

, _ . _ _ 12.45 3 , m 

leg P 14 4 /px 4 - - 3.8 - — — - - log T 8 

O 


(A25) 


(A26) 


52 


App. A - page 15 


This contribution to the total rate is dominant in the range 
0.4 < T < 1.7. At T > 1.7, the contribution from the 4.844 MeV 

o 8 

level becomes dominant. Using the same experimental data we 
. 2 

obtain u)0 = 0.05 and 

a 

22.18 

log P 14,4 /px 4 = 2 " “T 3/2 log T 8 (A27) 

1 8 

At temperatures Tg > 3, the rate obtained from table (\) should 
become valid. 


A- 4 0 16 (a. v)Ne 20 

Almqvist,E. and Kuehner, J.A. , and also Evans et al (to be 

published 1964) report the following results from measurements 

20 

of some properties of excited levels of Ne 
For the 5.63(3-) level 
Ty/T = 0.077 + 0.008 

r r 

r = 1.9 x 10~ 3 ev, (2-t+l) = (1.0 + 0.1) x 10" 3 ev 

r 
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For the 5.80 (1-) level 


Ty/T < 3 x 10 


-4 


The capture rate becomes 


. (2.1 + 0.2) x 10 To— 45 9/To 

P u „/pX_ = S-7* X 10 

16,4 4 

8 


"372" 


(A28) 


in the range 2 < T 0 < 8. 

O 


The photodisintregation rate 


= 10 12 * 15 x 10“ 28 * 4 / T 9 

-Lb , 4 4 


good 
.s wi 

valid to Tg < 4 (within 50%) . 


(A29) 


This last rate is within 10% at T <8 and is most likely 

8 


A- 5 O 18 (a, v)Ne 22 0 18 (g,n)Ne 21 

The rate can be computed from table (1) . The neutron 

branching ratio can be obtained by multiplying the rate of 
18 20 

O (a,y)Ne by the ratio of neutron to total emission ratio as 
discussed in Appendix B. 
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A- 6 Ne££ (a, v)Mq24 

The experimental data covers the range E c m = 2.073 to 
3.269 MeV, corresponding to E* = 11.387 to 12.583 MeV in the 
Mg^ compound nucleus in which 11 resonances are reported. 

The value of (y/D) was found to be 10^*^7 and s = 10^*^ 

Q 91 

as compared to S = 10 J when calculated with the optical model. 

A— 7 Ma 24 (a. v) Si 28 

Here we have information in the range E =1.31 MeV to 

^ c .m 

28 

3.73 corresponding to E* = 11.296 to 13.714 in the Si compound 
nucleus. We have § Q 2 (y/D) = 10^ 0 *^ and S = 10^®'^. The 
optical model S is 10^° * 5 ^. At temperature above T = 16 

O 

(where r < r ) we evaluate from the experimental data 

y (1 

(P /D) = C = 1.4 x 10" 5 . The best choice of E*, both from the 
Y 

experimental data and from the calculation over ihe sum of the 
resonances is E* = 1.63 MeV. 

This yield log P/px^ =4.15 - 82.1/T g - 1/2 log T g (A30) 

In Figure 15 the rate computed by summing over the reso- 
nances is shown, together with the low temperature approximation 
(from Table ( I )), and the high temperature approximation. As 
discussed before, the rate obtained by summing over the resonances 
is expected to underestimate the rate at both ends. In the graph 
the energy range where we have information on individual resonances 
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is shown with respect to the Gamow Peak. (Calculations by Jay Hauben) 


A- 8 C 12 + C 12 

The Chalk River group (Vogt, 1964) has recently studied very 
thoroughly (experimentally and theoretically) the nuclear re- 
actions between and a few light nuclei (C^, N 1 ^, O^ 6 ) . 

The elastic scattering and the capture cross section data have 
been improved to about 20% accuracy. The best fit to the data 
is obtained with the following set of parameters: V Q = 50 MeV, 

W Q = 10 MeV, R Q = 5 . 77f , a * 0.40f. In Figure (16) the capture 
cross section thus obtained is compared with the experiments. 
Good agreement is obtained if the resonances are averaged over. 

The lower range of the optical model calculation has been 
fitted by a formula of the form 


log o = 17.30 - log E - 


37.87 

-p72 * 


0.35 E 


(A31) 


The resulting value of a is plotted in Fig. (16) . The difference 
between the optical model values and eqn.(A31) is less than 3% 
in the range 3.5 to 5.0 MeV. In the same range the cross section 
varies by a factor of more than four hundreds. The thermonuclear 
reaction rate can be written as 


log i0 (P/pX c ) = 26.37 


36.55 (1+0.070 T 9 ) 
m 1/3 


1/3 


- 2/3 log T 


(A32) 


56 


App. A - page 19 


In the range of stellar interest (T^ ~ 0.7) the formula (A32) 
gives a result which is about a factor of two smaller than the 
result given by Reeves (1962) or by Fowler and Hoyle (1964) . 

The uncertainty on the result has been brought down by quite 
a large amount. In view of the fact that the experimental results 
are good to 20% accuracy; that these results are well reproduced 
by the optical model; that the eqn. (A3 3) matches within a few per- 
cent the behavior of the optical model curve up to 0.60 of the 
Coulomb barrier energy, it is difficult to see how the rate eqn. 
(A32) could be wrong by more than a factor of three. 


noon 

A- 9 Na (p.g)Ne We have detailed information for incident 

proton energy (center of mass) varying from 100 keV to 778 keV. 

The thermonuclear reaction rate has been computed by 

summing over the resonances and also by using the average 

reduced value of uur T /T as calculated from the first 14 

a p 

resonances. Using the method described in Appendix A we get 


our 


r e 2n? e gE = 10 2 * 60 


res 


u, P 


(A33) 


The Coulomb factor used in eqn. (A3 3) is the proton factor. We 
neglect the alpha Coulomb factor since the Q value for 
Na 23 (p,a)Ne 20 is already 2.3 MeV. This approximation should be 
reasonably good. The rate is then given by the formula 
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log p 23 ./p ^1 = 14 ‘ 15 “ 19 ' 43 T 8 1/3 " 2/3 lo 9 T 8 " °' 02 T q 

In the Figure 17 the lower curve represents the sum over 

individual resonances for the Na 22 (p,a) Mg 24 while the upper 

curve is computed from the formula. The sum over resonances 

represents accurately the thermonuclear reaction rate in the 

range T < 10. Above this temperature high resonances (for 
8 

which we have no information) start to contribute. Then the 
upper curve should be better. As we expect r ~ in this 
range, this curve probably overestimates the rate by a factor 
of two or so. 


(A34) 
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Production of neutrons through endothermic branching of 
thermonuclear reactions between charged particles. 

A. Fractional number of reactions talcing place with 
relative energies higher than a threshold energy E fc . 

r 00 p® 

F = j o (E ) n (E ) vdE/ j n (E ) a (E ) vdE 
Et o 

(B-l) 

= [J exp[-EAT , -a/E 1//2 ]dE]/[e _T AE q /tt/2] 

E t 

where a/E^ 2 = 2Trn, t = 3E Q /kT / , E 0 ^ 2 - akT'/2, AE Q = AEq/St , 

T' has been defined in Appendix A. Here the integral in the 
denominator has been evaluated with the usual method of re- 
placing the integrand by a properly adjusted Gaussian. The 
treatment of the numerator requires a little more attention. 

We now define E/E Q = x, E^/E q = x t = 1 + u^ (we consider only 
u t >o) , and y = x fc - x. Then we develop the integrand in the 
numerator of eqn. (1) in a series in power of y (including 
terms up to y^) 

eat + a A 1/2 - x + -Itj + y [i- + y 2 | + ... (b-2) 

t t ' x t 
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Now inserting this development in eqn. (1) and integrating 


we get, 

5/4 

x t 

F = — [1 - cp(z t ) 1 exp 
2 ^ 



(B-3) 


where 


Z t = / t/ 3 (x fc 5/4 



) 


and 



(2//tt) exp (-z 2 )dZ 


as tabulated e.g., in Jairike-Emde , page 24. 


Useful approximations : 


for u.«l, (x 5/2 _ 5x — + 9) ~ (15/8 u 3 _ 45/32 u 

t t t yx^ — t t 


and 


F ~ [1 - cp(Z t ) ] 7 


*kT 

for u ^ 0.5, F ~ — 
t AE, 


expf-^(x t + 2//x fc - 3) J 


(1 - X 


(B-4) 


B. Average number of neutrons emitted per collision. 
Here we need 


g = J°° a (E)n(E)vdE (E)n(E)vdE . 


(B-5) 


where r n /r is the ratio of the neutron width r to the total width 
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of the levels of compound nucleus, and E , here, becomes the 

t 

threshold energy for neutron emission. The ratio (r /T) is 
small only for (E-E^_) less than a few keV. In the rest of the 
energy range its value remains roughly constant. Depending upon 
the reaction under consideration it takes values ranging from 
(r n /D 1 to ~ 0.1. Neglecting the lowest part of the range 
the average number of neutrons emitted per collision becomes 

g ~ (iV r) F 

Ip 21 

0 (g,n)Ne Q = -700 keV. In the range of interest (r /T) ~ 1 

since r is still fairly small because of Coulomb effects and 
a 

no other channel is open (Q(a,p)~ -5.7 MeV) . As a function of 
temperature we obtain 


T 

8 

X 

t 

Z 

t 

F 

g 

o 

• 

CN 

1.78 

3.27 

~10 -3 

~10 -3 

3.0 

1.36 

1.38 

0.05 

0.05 

3.5 

1.22 

0.85 

0.17 

0.17 

4.0 

1.12 

0.44 

0.31 

0.31 

Here g is also 

i p 

the number of neutrons produced per each O 

nucleus in the 

gas . We 

notice 

that at 

higher temperature this 


21 24 

number will be doubled by the Ne (a,n)Mg reaction. 
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22 , 25 

Ne ( g,n)Mg 


Q = -482 keV 


Par the same reasons as in the previous case we expect 
(r n /D ~ 1 as we have Q(a,p) = -3.5 MeV. 


T 

8 

X 

t 

Z 

t 

F 

g 

1.5 

1.26 

1.23 

0.07 

0.07 

ro 

• 

o 

1.04 

0.20 

0.38 

0.38 

the 

number of 

neutrons 

per Ne 22 

nuclei. Higher 


temperature may double the crop of neutrons if enough He are 
still present by then. 


C 12 (C 12 , n)Mg 23 Q =*2.603 MeV. 

The experimental evidence (Bromley 1960) gives a branching 
ratio of 0.05 (fairly constant) from 10 MeV down to 5 MeV in the 
center of mass system. 

In the lower energy range two channels are still widely 
open r (C l2 ,a) ; Q = 4.619 MeV, (C 12 ,p); Q = 2.230 Mev] . It is 
probably reasonable to assume that the value (r n /r) = 0.05 is 

valid at lower energies, although this may be uncertain by a 
factor of two or three. 

To evaluate F we use E = 2 .418 (T' ) 2//3 , where T' = 

o 9 9 

Tg/(1+ 0.070 Tg) , as discussed in Appendix A. 
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T 

X 

Z 

F g/2 

F 

% 

9 

t 

t 

app 

exact 

* 

0.6 

1.56 

2.82 

5.9xl0” 4 2xl0" 5 

5 . 5xl0 -4 


0.7 

1.41 

2.02 

8JSxlO“ 3 1 . 5xl0~ 4 

7. lxlO -3 


0.8 1.30 1.41 0.045 lxlO" 3 0.036 

0.85 1.25 1.18 0.079 2xl0 -3 0.065 

0.9 1.20 0.97 0.12 3xlO -3 0.105 

1.00 1.13 0.58 0.24 0.005 0.215 

1.10 1.06 0.28 0.37 0.01 0.373 

1.20 0.012 0.487 

As each collision tc&es two C , we have listed in the last 

column g/2, the number of neutrons per in the gas. In the 

last column the exact value of F (computer integrated) is given 

(calculations by E. Milford). 
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Just as an illustration we may consider five red giant 
stars which show very different composition. In the table the 
ratio of iron to hydrogen (N^/N^, related to the activity of 
the e process) ; the ratios of Ba-La-Ce-Nd to Fe, (N/Np e , related 
to the activity of the s process) and the ratios of C and 0 to 
Fe, (possibly related to the moment of occurrence of the 
s process) are given. The solar abundances are used as a 
standard. 


Star 

(N Fe /N H> 

(N /N ) 
s' Fe 

(N /N ) 
cr Fe 

(N /N ) 
o' Fe 

Sun 

1 

1 

1 

1 

HD122563 

~ 0.01 

0.02 



HD165195 

~ 0.01 

~ 1 



HD201626 

0.03 

15 

5 

< 1.6 

HD26 

0.2 

15 

5 

< 0.5 

e Virginis 

1 

1 

1 

1 


References: Wallerstein (1964), Wallerstein et al (1963). 


The question of whether evolved material can be brought to the 
surface of a star during its quiet lifetime is by no means 
settled. In this discussion we shall naively assume that such 
transport takes place at least in red giants. This assumption 
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is not very important. A similar discussion could be made on the 
content of the interstellar gas from which these stars were made 
if the transport assumption is wrong {or crude) . 

The ratio of N /N is thought to be related to the age of 
iTG H 

the star. The two first stars (apparently very old) are quite 
similar in many abundances (not quoted here) except for their s 
process abundance. This may reflect the fact that as shown 
above, neutron processes are not an automatic consequence of 
helium stage (and of more advanced stages) but depend amongst 
other things upon the stellar masses. The two following stars 
present evidence of intense neutron activity (increase in the 
s process elements) . The enhancement of C without enhancement O 
suggests that the process took place in the early days of helium 
burning, possibly from proton admixture or from burning at 
high temperature in the peak of the helium flash. 

Finally e Virginis is given as an example of ’'normal" red 
giants, in many ways similar to the sun. 

Let it be emphasized that these analyses (and others in the 
text) are very preliminary and serve mostly as illustrations. 
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Table of Q Values for Various Reactions 


Target 


Q Values MeV 




(P/ Y) 

(p,°0 

(P,n) 

(a, y) 

(a, n) 

c 12 

1.943 

-7.557 

-18.390 

7.161 

-8.507 

c 13 

7.549 

-4.064 

-3.004 

6.357 

2.214 

n 14 

7.291 

-2.922 

-5.931 

4.404 

-4.737 

N 15 

12.125 

4.964 

-3.543 

4.012 

-6.430 

0 16 

0.598 

-5.218 

— 

4.730 

-12.145 

017 

5.597 

1.193 

-3.544 

7.347 

0.588 

0 18 

7.992 

3.980 

-2.450 

9.667 

-0.700 

f 19 

12 . 844 

8.114 

-4.031 

10.465 

-1.949 

Ne 20 

2.454 

-4.132 

-16.100 

9.314 

- 7.221 

Ne 21 

6.743 

-1.750 

-4.305 

9.885 

2.555 

Ne 22 

8.790 

-1.675 

-3.624 

10.616 

-0.482 

Na 23 

11.693 

2.379 

-4.842 

10.098 

-2.971 

Mg oc 

2.287 

-6.859 

-14.800 

9.986 

-7.193 

Mg 25 

6.301 

-3.142 

-5.044 

11.133 

2.655 


8.272 

-1.826 

- 4.797 

10.650 

0.036 

Al 26 

7.471 

-1.871 

— 

10.417 

-0.910 

Al 27 

11.581 

1.594 

-5.598 

9.664 

-2.652 

Si 2 ® 

2.735 

-7.700 

-14.580 

6.946 

-8.135 

Si ,n 

5.585 

-4.832 

-5.743 

7.111 

-1.532 

Si 30 

7.286 

-2.378 

-5.030 

7.917 

-3.504 

p3 1 

8.863 

1.917 

-6.218 

6.998 

-5.572 

s 32 

2.285 

-4.211 

-13.780 

6.641 

-8.628 

s 33 

5.218 

-1.526 

-6.358 

6.792 

-2.002 

c 34 

& o tr 

6.367 

-0.631 

-6.203 

7.213 

-4.629 

cl 35 

8.506 

1.865 

-6.763 

7.213 

-5.866 

s 36 

8.402 

0.544 

-1.920 

6.809 

-3.066 

A 36 

1.890 

-4.355 

— 

7.044 

-8.685 
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Figure 1 . Iso-intensity curves for emission of neutrinos in 
the p-T plane. The label on each curve identifies the 
energy yield in erg/gm/sec along that line. Three processes 
are considered: plasma neutrinos, photoneutrinos, pair 

annihilation neutrinos. The dashed parts of the curves 
reflect the fact that photoneutrino rates in degenerate 
matter have not yet been properly computed. 
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Figure 2 . Densities and temperatures in hydrogen-burning core 
and shells for stars of 1.09, 4 and 15.6 solar masses The 
numbers on the curves identify the length of time elapsed 
since the beginning of the hydrogen-burning phase and, 
further out on the curves the length of time since the 
onset of the hydrogen burning shell. The results are from 
Hayashi et al (1962) . 
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Figure 2. 





Figure 3 . Rates of various nuclear reactions of importance 
during the helium burning phase. The electron screening 
factors are not included here. This effect will be 
important mostly during the period preceding the nitrogen- 
helium flash in small stars. 
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figure 4 . Behaviour of various quantities during the pre- 
helium flash period of a 1.3^M 0 star. The calculations 
have been made using a sequence of model by Schwarzschild 
et al (1962) . L/M is the mean energy generation rate for 
the whole star (mostly from the hydrogen burning shell) . 
L core/ M core (gravitation) is the energy generation rate from 
the contraction of the helium core, L v (core) is the rate of 
neutrino dissipation by plasma neutrino processes, core 

is the nuclear energy released by the 3He 4 -* C ^ burning. 

€^ 2 ' the energy from the cl3 (He 4 , n) reaction assuming 
X]o = 10“ 4 (Pop I star), "e '^4 the energy from the N^ 4 (He 4 , y)F^® 

( 8 , vJoiS assuming X 14 = 10~2 (Pop I star) or X 14 = 10 “ ^ 

(Pop II star) . The effect of these two last nuclear reactions 
on the sequence of models have been neglected. (Calculation 
by A. Liebman) 
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Figure 



MODEL NO ( SCHWARZSCHILD) 



Figure 5 . Energy dissipation rate from plasma neutrinos in 
erg/gm/sec as a function of T g (as a parameter) and 
log p (as the ascissa) . 
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Figure 6 . Densities and temperatures in helium-burning core 
and shells arid carbon-burning core for stars of 0.7, 4 and 
15.6 M 0 . The numbers on the curves identify the fractional 
weight of helium left in the core. The calculations pertain- 
ing to the carbon-burning phase do not take into account the 
effect of neutrino dissipation processes. The results are 
from the work of Hayashi et al (1962) . 
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Figure 




Figure 7 . Iso-abundance curves describing isotopic outcome of 
helium burning processes as a function of (fixed) density 
and temperature. Except for Neon, only the most important 
species is mentioned at any point of the p-T plane. The 
caption e,g, c 12 >50% in a given region means that in the 
physical conditions pertaining to that region the 
fractional weight of C^ 2 is more than 50% (but less than 80%) . 
The rates used in this figure differed slightly from the 
rates quoted in the text. Present rates would allow Ne 2 ® 
to go to about 20% in the region where Ne 2 ® ^ 10% is quoted. 
Everywhere else the abundance of Ne 2 ® is at best a few 
percent . 
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Figure 
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Figure 8 . Iso--abundance curves describing the isotopic outcome 
of the helium burning phase as a function of stellar masses. 
(Salpeter and Deinzer (1964) ypure helium stars) . Three 
curves are given for C 12 and 0 16 , representing the choice 
of 0.1, 0.4 and 1.0 for the value of the parameter 0^ 

(see text) . The values of Ne 2 ® and Mg 2 ^ are practically 
independent of the choice of 9^ . 
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Figure 





Figure 9 . Energy dissipation by pair annihilation neutrinos 
in erg/gr/sec. The alscissa is and the parameter log p 
en gm/ cm^ . 
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Figure 9. 
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Figure 10 . Carbon-burning stage curves. The "first model" 
line describes the central stellar conditions where the 
nuclear energy generation rate (from C ^2 + ) equals 

the neutrino energy dissipation rate in the center of the 
star. The "corrected model" curve described the central 
stellar conditions when the total (volume integrated) 
nuclear energy generation rate equals the total neutrino 

energy dissipation rate (assuming an n = 3 'polytropic 
model) . 
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Figure 11 . Lifetime and energy generation rate of the carbon 
burning stage as a function of the central temperature 
using the values given in Figure 10 and described in the 
text. 
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Figure 11 . 




Figure 12 . Oxygen and neon -burning stage curves. The model 
is described in the text and in the caption of Figure 10. 


91 



Figure 12. 
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Figure 13 . Graph of £ (2-1+1) g^ /g 0 as a function of y ' 

(see eqn. A3). The graph is valid for y 1//2 > 4, and when 
the relative energy is less than one-half of the Coulomb 
barrier energy (B = Zj^e^/R) ♦ 
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Figure 14 . Values of the parameters S and g (see App. A) as a 
function of the target nucleus for (a, y) reactions The 
points refer to the results of optical model calculation, 
to the matching formula eqns. (A5) and (A7) , and for Ne^®(a,y) 
and Mg^(a,y) to a determination of the strenath function 
from the level scheme of the compound nucleus. (Calculations 
by A. Weiswasser) 
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Figure 15 . Rate of Mg^^(a,y)Si^ 

The sum over individual resonances ( ) is expected to be 

best in the range from 6 < Tg < 10. Below and above this 
range unanalyzed resonances are expected to play a role 

hence to increase the rate. The curve called optical model ( ) 

is best in the range Tg = 3 to about 10. Below this temperature 
the randomness in the position of the levels makes of it an 
upper limit. Above this temperature the fact that T a becomes 
larger than T makes again of it an upper limit, (the optical 
model yields only the total capture rate). The intermediate 
curve in the higher temperature range is based on the experi- 
mental gamma strength function. It should be best there. 






1 o 12 

Figure 16 . The C + C total reaction rate. 

The ( — - ) curve represents the experimental results. 
The solid curve is the result of the optical model. The 
dashed curve (merging in the solid curve at low energy) is 
obtained from the fitting equation no 
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Figure 17 . Rate of the Na^(p,a)Ne 

The dashed line ( ) is the sum of individual rates. 

It should be correct in the range Tg < 10. Missed 
resonances make of it a lower limit at higher temperature. 
The dotted line ( . . ) has been obtained with the equi- 
valent of a strength function. It overestimates the rate 

at Tg > 15 . There a best choice rate is quoted ( ) . 

(Calculation by Jay Hauben) 


log P/px i 




ACKNOWLEDGMENT 


I acknowledge the help of Messrs. J. Hauben, 
A. Liebman and E. Milford, and Miss A. Weiswasser 
for numerical calculations; Mr. T. Psaropulous in 
drawing the figures; and Mrs. Y. Lecavalier and 
M. Yastishak in typing the manuscript. 

I wish to thank Dr. R. Jastrow and the staff 
at the Institute for Space Studies for their 
hospitality. 


